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Об одном подходе к реализации 
автоматического контрастирования 
тепловизионного изображения 


В статье рассматривается один из способов реализации автоматического контрастирования тепловизионного 
изображения в специализированных вычислительных системах, предназначенных для предварительной 
обработки сигналов матричных фотоприемных устройств. Предлагаемый способ автоконтрастирования 
основан на сочетании метода линейного преобразования диапазона входных кодов пикселей в диапазон 
выходных кодов и метода эквализации гистограмм и ориентирован на аппаратно-программную реализацию. 


Введение 


Матричные фотоприемные устройства (МФПУ) широко применяются в качестве 
приемников инфракрасного излучения различных тепловизионных систем технического 
зрения, обеспечивая решение широкого круга задач оборонного и общегражданского 
назначения. Разработка математических методов, новых алгоритмов и аппаратно - 
программных средств, предназначенных для обработки информации, получаемой от 
этих фотоприемников, является важнейшей задачей, характерной для периода разви- 
тия современных поколений приборов инфракрасной фотоэлектроники [1]. 

Немаловажную роль в процессе обработки и визуализации тепловизионной инфор- 
мации играют специализированные вычислительные системы, выполняющие преобра- 
зование аналоговых сигналов от фоточувствительных элементов МФПУ в цифровую 
форму, формирование тепловизионного изображения на основе полученной информа- 
ции и его предварительную обработку, вывод этого изображения на экран черно-белого 
монитора. 

Одним из основных алгоритмов предварительной обработки тепловизионного 
изображения является алгоритм автоматического контрастирования изображения. 


Постановка задачи 


Современная технология изготовления МФПУ не позволяет без предваритель- 
ной обработки получить тепловизионное изображение, пригодное для визуального 
восприятия. Основная проблема заключается в существовании геометрического шума 
и неоднородности изображения, связанных с различной чувствительностью и темно- 
выми токами отдельных фоточувствительных элементов (ФЧЭ) матрицы. В связи с 
этим при обработке сигналов ФЧЭ возникает необходимость их коррекции по зара- 
нее проведенным калибровкам. Типовая коррекция проводится по формуле 
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где $,, — скорректированный элемент изображения от 1, -го ФЧЭ матрицы; 
5, ; — входной, некорректированный элемент изображения от 1, }-го ФЧЭ матрицы; 
КТ, ‚ — значение калибровки по темному фону для 1, -го ФЧЭ матрицы; 


К2, ‚ — значение калибровки по светлому фону для 1, -го ФЧЭ матрицы. 


Очевидно, что подобная коррекция может быть достаточно просто проведена в 
цифровом виде. Для различных матриц и разных калибровочных значений измене- 


ние 5, может быть весьма значительным, в то время как интересующая оператора 


динамика яркостей по сцене, воспринимаемая человеческим глазом, намного меньше. 
В связи с этим возникает необходимость выделения из всего возможного диапазона 
кодов, получаемых после применения коррекции, только того диапазона, в котором 
сосредоточено текущее распределение яркости по наблюдаемой сцене и преобразо- 
вания этого диапазона к фиксированному, воспринимаемому человеческим глазом — 
8 бит, или 256 градаций. 

Подобная задача автоконтрастирования изображения может решаться несколь- 
кими способами. Самый простой из них заключается в выделении 8 старших 
разрядов после применения коррекции. Этот способ применим только в том случае, 
когда калибровки КП и К2 имеют строго определенное фиксированное значение 
вблизи соответственно минимума и максимума возможных яркостей и динамика 
яркостей по сцене всегда перекрывает всю шкалу возможных значений 5 из 
выражения (1). В любом другом случае подобный метод приведет если не к полной, 
то к частичной потере полезного сигнала. Так, если изображение с матрицы ФЧЭ имеет 
динамику 48 дБ(8 разрядов) или менее, а после коррекции диапазон кодов, например, 
составляет 96 дБ (16 разрядов), то простое выделение 8 старших разрядов полностью 
уничтожит информацию об изображении, оставив только среднюю по кадру яркость. 

Следующим способом (по сложности обработки) является линейный пересчет 
диапазона входных кодов в диапазон выходных кодов по пиковым значениям яркости в 
кадре [2]. В этом случае находятся динамические уровни черного и белого как 
минимальное и максимальное значения сигнала от всех ФЧЭ матрицы за один кадр. 
Разность этих значений и определяет масштабный коэффициент, применяемый для пе- 
ресчета входных значений в выходные. Однако применение этого метода на практике 
ограничено некоторыми условиями. Как известно, матрицы ФЧЭ могуг содержать 
дефектные элементы, сигналы от которых не соответствуют их освещенности и могут 
приниматься за ложные уровни черного и белого. Кроме того, свою роль при слабой 
освещенности могут сыграть и собственные шумы ФЧЭ, и короткие помехи, наведен- 
ные в аналоговом канале и отображающиеся на выходном изображении как одиночные 
пиксели, отличающиеся по яркости. Все это приводит к тому, что канал преобразования 
данных от матрицы до видеомонитора «слепнет», реагируя на неверные уровни черного 
и белого. Существует также и объективная причина «ослепления» канала — очень яркие 
(темные) объекты, попадающие в поле зрения матрицы, но не представляющие 
интереса для наблюдателя. Например, попавшая в поле зрения матрицы спираль лампы 
накаливания сделает неразличимо темными окружающие слабоизлучающие предметы. 

Третий способ, называемый методом эквализации гистограмм [2], является ста- 
тистическим. Метод заключается в разбиении всего диапазона кодов после коррекции 
(исходного диапазона кодов) на участки по количеству выходных градаций яркости. 
Причем разбиение производится таким образом, чтобы все участки содержали оди- 
наковое количество пикселей кадра. Очевидно, что в общем случае кодовая длина 
каждого диапазона не будет равна остальным. Далее всем пикселям, принадлежащим 
определенному участку, ставится в соответствие единственное значение кода выход- 


544 «Искусственный интеллект» 372009 


Об одном подходе к реализации автоматического контрастирования... 


10к 


ного изображения. На практике этот метод применим к статическим или медленно 
меняющимся изображениям. Обработка же быстро меняющихся изображений требует 
огромного объема вычислений в реальном времени (в случае программной реализа- 
ции алгоритма) или больших аппаратурных затрат (при аппаратной реализации). 

Таким образом, разработка алгоритма автоматического контрастирования теп- 
лового изображения, обеспечивающего качественную визуализацию тепловых объектов 
и одновременно не предъявляющего высоких требований к аппаратным средствам и 
вычислительной мощности специализированных вычислительных систем, является весь- 
ма актуальной. 


Реализация автоконтрастирования тепловизионного 
изображения, ориентированная на аппаратно- 
программную реализацию 


Предлагаемый вариант реализации автоконтрастирования изображения сочетает 
в себе относительную простоту первых двух методов автоконтрастирования и преи- 
мущества статистического подхода, позволяющие избежать случайного «ослепления» 
канала видеоданных. Алгоритм основан на применении не пиковых, а интегральных 
динамических уровней черного и белого в сочетании с переключением поддиапазо- 
нов используемых разрядов значений пикселей. На рис. 1 представлена структура 
аппаратных средств, реализующих алгоритм автоконтрастирования. 


$ п [7 > 
от схемы № 
‘оррекции $ БУ 
=> 
БСЯ к — 
ы 
данные 
> тистограммы й 
„„  БФГ "о МК БФПИ » 
7 нЕ ле [АП 


Рисунок 1 — Структура аппаратных средств, реализующих алгоритм 
автоконтрастирования тепловизионного изображения 


На рис. | представлены следующие основные блоки (будут рассмотрены под- 
робнее ниже): 

БСЯ - блок определения средней по кадру яркости; 

У — вычислитель модуля яркости пикселя, производит также вычитание средней 
по кадру яркости из значения яркости текущего пикселя; 

Мх - мультиплексор, производящий предварительный выбор значащих разрядов 
яркости пикселей; 

УУ - узел управления мультиплексором; 

БФГ - блок формирования гистограммы; 

МК - микроконтроллер, производящий обработку гистограммы и вычисление ло- 
кального коэффициента умножения; 
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БУ - блок умножителя; 

К - компандер; 

БФП - блок формирования пикселя выходного изображения. 

В качестве динамических уровней черного и белого в рассматриваемом алго- 
ритме автоконтрастирования применяются значения гистограммы распределения яркос- 
ти по кадру, ниже и выше которых лежит заданное количество М пикселей в кадре. 
Это позволяет полностью исключить влияние на динамические уровни черного и бе- 
лого такого яркого (темного) объекта, размер которого на изображении в пикселях 
меньше М. 

Гистограмма распределения яркости по кадру и, соответственно, динамические 
уровни черного и белого вычисляются максимально точно, когда количество поддиа- 
пазонов гистограммы (точек по оси абсцисс) равно количеству кодов входных пикселей 
после коррекции. При этом, например, если результат коррекции имеет 12 разрядов, 
необходимо организовать 4096 поддиапазонов, что приведет к огромным затратам 
оборудования при аппаратной реализации формирования гистограммы. Поэтому для 
аппаратной реализации алгоритма при вычислении гистограммы распределения яркостей 
пикселей весь диапазон входных кодов был разбит на 64 поддиапазона с использова- 
нием дополнительной регулировки переключением разрядов, что является оптимальным 
с точки зрения аппаратных затрат и адекватности регулировки контраста изображения. 

Для переключения используемых разрядов и вычисления гистограммы приме- 
няется не исходное значение яркости пикселя после коррекции, а его модуль (далее — 
модуль яркости), что позволяет работать также с пикселями, яркость которых ниже 
значения калибровки по темному фону К1. Узел управления мультиплексором УУ 
анализирует наличие пикселей, модуль яркости которых имеет старшие значащие раз- 
ряды МЪЗВ, (МЗВ-1), (МЗВ-2), (МЗВ-3), а также количество таких пикселей в кадре. 
На основании превышения (непревышения) определенного порога (например, 1% от 
общего количества пикселей в кадре) принимается решение о выборе одного из 4 диа- 
пазонов кодов модуля яркости: [0 — (М$ЗВ-3)], [1 — (МЗВ-2)]|, [2 — (М$В-1], [3 —- МЗВ]. 

Выбранный диапазон кодов разрядности (п-2) после мультиплексора разбит на 
64 поддиапазона. В процессе обработки данных кадра изображения, вводимого с мат- 
рицы ФЧЭ, происходит накопление информации о количестве пикселей, модули яр- 
кости которых попадают в соответствующие поддиапазоны блоком формирования 
гистограммы БФГ. Далее последовательным суммированием количества пикселей по 
поддиапазонам «снизу вверх» и «сверху вниз» в микроконтроллере МК вычисляются 
динамические уровни черного и белого — как номера поддиапазонов, ниже и выше 
которых соответственно находится не менее М пикселей кадра. На основе полученных 
динамических уровней черного и белого микроконтроллер вычисляет локальный коэф- 
фициент умножения К,, который может принимать значения от 1 (если уровень чер- 
ного равен 0 — первый поддиапазон, а уровень белого — 63 — последний поддиапазон) 
до 64 (если уровень черного равен уровню белого, т.е. весь диапазон яркостей сосре- 
доточен в одном поддиапазоне гистограммы). В общем виде значение коэффициента 
умножения можно выразить следующим образом: 

К, ) 
тах— шш-+1 
где тах — уровень белого, тт — уровень черного. 

Необходимо отметить, что для правильного применения локального коэффициен- 
та умножения необходимо, чтобы среднее за кадр значение яркости пикселей, умно- 
жаемых на этот коэффициент, равнялось нулю. Это обеспечивается применением 
аппаратного вычислителя средней по кадру яркости и сумматора-вычитателя »', вы- 
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читающего среднее значение из яркости каждого пикселя. После умножения модулей 
яркости пикселей на локальный коэффициент умножения яркости всех (за исключе- 
нием 2М) пикселей, будут лежать в полном диапазоне выходных кодов - 7 бит. Для 
сохранения информации о яркости оговоренных 2М пикселей в диапазоне выходных 
кодов предусмотрены защитные интервалы сверху и снизу (для модуля яркости этот 
защитный интервал необходим только сверху) и компандирование: яркость основной 
массы пикселей линейно трансформируется в диапазон между нулем и защитным ин- 
тервалом, а 2М пикселей, не попадающие в основную массу, — с загибом амплитудной 
характеристики в диапазон от верхнего защитного интервала до максимума шкалы 
выходных кодов. После компандера в блоке формирования пикселя выходного изоб- 
ражения БФП производится восстановление исходного знака яркости пикселя и до- 
бавление желаемой средней яркости, равной половине шкалы выходных кодов, что 
повышает субъективную различимость изображения. 


Заключение 


Предлагаемый способ автоконтрастирования был реализован аппаратно-программ- 
ным способом в многопроцессорной вычислительной системе, предназначенной для 
обработки сигналов МФПУ [3]. Вся схема, представленная на рис. 1, кроме микроконт- 
роллера, реализована аппаратными средствами ПЛИС, а анализ гистограммы, нахожде- 
ние динамических уровней черного и белого, вычисление локального коэффициента 
умножения выполняется программным способом. В вычислительной системе нормаль- 
ный (единичный) коэффициент усиления канала видеоданных соответствует ситуации, 
когда выбран второй вход мультиплексора [1...(п-2)] и локальный коэффициент уси- 
ления равен единице. Таким образом, имеется возможность уменьшить сигнал до 
восьми раз (ситуация, когда калибровочные значения меньше полезного видеосигнала) 
или усилить до 128 раз (полезный видеосигнал много меньше калибровочных значе- 
ний). Общая динамика коэффициента усиления видеотракта составляет 60 дБ. 


Литература 


1. Пономаренко В.П. Инфракрасная техника и электронная оптика. Становление научных направ- 
лений (1946 — 2006) / В.П. Пономаренко, А.М. Филачев. - М. : Физматкнига, 2006. - 336 с. 

2. Гонсалес Р. Цифровая обработка изображений / Р. Гонсалес, Р. Вудс. — М. : Техносфера, 2005. - 1072 с. 

3. Катаев О.В. Многопроцессорная вычислительная система обработки сигналов матричных фото- 
приемных устройств / О.В. Катаев [и др.] // Многопроцессорные вычислительные и управляющие 
системы -— 2007 : материалы Международной научно-технической конференции. Т. 1. — Таганрог : 
ТТИ ЮФУ, 2007. — С. 113-118. 


О.В. Катаев, П.А. Соломахи, А.В. Бовкун 

Про один шдх до реалзацй автоматичного контрастування теплов1з ного зображення 

У статти розглядаеться один 13 способв реалзацй автоматичного контрастування теплов1з1йного 
зображення у спещалзованих обчислювальних системах, призначених для попередньо! обробки сигнал!в 
матричних фотоприймальних пристро!в. Пропонований спос1б автоконтрастування заснований на 
поеднанн! методу лийного перетворення д1апазону вхадних код1в шкселв у дапазон вихдних код\в 1 
методу еквал1зацй г1стограм 1 ортентований на апаратно-програмну реализацию. 


О.Т. Кашуеу, Р.А. ЗоотаЕйт, А.Т. ВоуКип 

АБощ Фе Арргоасв 0 ВеаН7айоп о? 'ТВегта|! Ппасе Ашотайс Зато 

ш Фе рарег 15 415си5зе4 опе те@фо@ оЁ Феппа| ппазе ащютайс зашше ш зрес1а]-ригрозе сотршег 
зу$бетз УСН аге ииепде4 Юг $1епа[$ ргергосеззт» ш таблх рВоюдаесюгз. ТБе зисоезе4 тейо4 оЁ 
ащютайс зашише 15$ Базед оп сот таНоп оЁ мо тео4д$. ТБе Йг$ опе 15 тефо4 оЁ Ппеаг тапзЮгтаНоп 
оЁ шриё р!хе|5 содез гапее шо ошёрий соде$ гапсе. ТВе зесоп4 опе 1$ те@фо4 оЁ Б15юэгатз едиаНтайоп. 
ТБе тефоа оЁ ащютайс запиио 1$ опете4 оп Ваг4ууаге ап4 зоЙ\уаге геай тай оп. 
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«Штучний 1нтелект» 372009 547 


